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Abstract 

Correlation effects on the electronic structure of half-metallic CrAs in 
zinc-blende structure are studied for different substrate lattice constants. 
Depending on the substrate the spectral weight of the non-quasiparticle 
states might be tuned from a well developed value in the case of InAs 
substrate to an almost negligible contribution for the GaAs one. A piezo- 
electric material that would allow the change in the substrate lattice pa- 
rameters opens the possibility for practical investigations of the switch- 
able (tunable) non-quasiparticle states. Since the latter are important for 
the tunneling magnetoresistance and related phenomena it creates new 
opportunities in spintronics. 
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I. INTRODUCTION 



One of the strongest motivations to investigate magnetic semiconductors and half- 
metallic ferromagnets (HMF) is the possibility to design and produce novel stable 
structures on semiconducting substrates with new interesting properties. Adopting 
this point of view first-principle studies are an excellent starting point to predict new 
systems having the desired properties. Recently Akinaga et al. [1] found the possibility 
to fabricate zinc-blende (ZB) type CrAs half-metallic ferromagnetic material. Experi- 
mental data confirmed that this material is ferromagnetic with the magnetic moment 
of 3/XB) in agreement with theoretical predictions [1]. According to this calculation this 
half-metallic material has a gap of about 1.8eV in the minority spin channel which 
has attracted much attention to this potential candidate for spintronic applications, 
keeping in mind also its high Curie temperature T c around 400JT. Note that recent 
experiments on CrAs epilayers grown on GaAs(OOl) evidenced an orthorhombic struc- 
ture, different from the ZB one, so the structure is rather sensitive to the preparation 
process [2]. However, it is highly desirable to explore the possibility of existence of 
half-metallic ferromagnetism in materials which are compatible with important III-V 
and II-IV semiconductors. For this purpose efforts has been made on the metastable 
ZB structures, such as CrAs [1,3]. Beyond the importance of the preparation tech- 
niques, attention should be given to the understanding of the finite temperature prop- 
erties in these HMF materials. Therefore, it is interesting to explore theoretically the 
mechanism behind half-metallic ferromagnetism at finite temperature from a realistic 
electronic structure point of view. 

Theoretical studies [4] of the 3d transition metal monoarsenides have shown that 
the ferromagnetic phase of ZB structure CrAs compound should be more stable than 
the antiferromagnetic one. The calculated equilibrium lattice constant is larger than 
that of GaAs (5.65A) having a value of 5. 80 A [4]. Following this work similar electronic 
structure calculations concerning the stability of the half-metallic ferromagnetic state 
in the ZB structure have been carried out [5]. 

However standard local density approximation LDA (or GGA, LDA+U, etc.) calcu- 
lations are in general insufficient to describe some important many-particle features of 
the half-metallic ferromagnets. One of these many-body features, the non-quasiparticle 
(NQP) states [6-8] contribute essentially to the tunneling transport in heterostructures 
containing HMF [9,10]. The origin of these states is connected with the "spin-polaron" 
processes: the spin-down low-energy electron excitations, which are forbidden for HMF 
in the one-particle picture, turn out to be possible as superpositions of spin-up elec- 
tron excitations and virtual magnons [6,7]. Recently we have applied the LDA+DMFT 
(dynamical mean field theory) approach [11,13] to describe from first principles the 
non-quasiparticle states in a prototype half-metallic ferromagnet NiMnSb [14]. 

In this paper, we describe the correlation effects in CrAs HMF material in the 
framework of the LDA+DMFT approach. We will show that these many-body spin- 
polaron processes are very sensitive to structural properties of the artificially produced 
CrAs compound. Depending on the substrate characteristics, such as a large lattice 
constant, as in the case of InAs (6.06A), or a smaller one as in the case of GaAs 
(5.65A) the spectral weight of the non-quasiparticle states can be tuned from a large 
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value in the former case (InAs) to an almost negligible contribution in the later case 
(GaAs). Therefore, the correlation effects in conjunction with structural properties 
determine the behavior of electronic states in CrAs near the Fermi level which has a 
substantial impact on the tunneling transport in the corresponding heterostructures. 

II. TUNNELING TRANSPORT AND NON-QUASIPARTICLE PHYSICS 

First, let us explain in a simple way why the non-quasiparticle states are important 
for the tunneling transport (for formal derivations, see Refs. [9,10]). To this aim we 
consider the case of a narrow-band saturated Hubbard ferromagnet where the current 
carriers are the holes in the lowest Hubbard band and the non-quasiparticle states pro- 
vide all spectral weight for the minority spin projection [7]. A schematic density of 
states is shown in Fig. la. Suppose we have a tunnel junction with two pieces of this 
ferromagnet with either parallel (Fig. lb) or antiparallel (Fig. lc) magnetization direc- 
tions. From the one-particle point of view, the spin-conserving tunneling is forbidden in 
the latter case. However, in the framework of many-particle picture the charge current 
is a transfer process between an empty site and a single-occupied site rather than the 
motion of the electron irrespective to the site like in the band theory and therefore the 
distinction between these two cases (see Fig. 1), is due only to the difference in the den- 
sities of states. It means that the estimations of the tunneling magnetoresistance based 
on simple one-electron picture is too optimistic; even for antiparallel spin orientation 
of two pieces of the half-metallic ferromagnets in the junction for zero temperature the 
current is not zero, due to the non-quasiparticle states. More exactly, it vanishes for 
zero bias since the density of NQP states at the Fermi energy equals to zero. However, 
it grows with the bias sharply, having the scale of order of typical magnon energies, 
that is, millivolts. 

The latter statement is confirmed by a formal consideration of the antiparallel spin 
orientation case, based on the standard tunneling Hamiltonian [15]: 

H = H L + H R + E( T fe P 4 T c pi + h.c.) (1) 

kp 

where H L R are the Hamiltonians of the left (right) half-spaces, k and p are correspond- 
ing quasimomenta and spin projections are defined with respect to the magnetization 
direction of a given half-space with opposite magnetization directions (spin is supposed 
to be conserving in the "global" coordinate system). For the tunneling current I in the 
second order in T kp one has 

Joc]T|T fep | 2 [l + A g -/(V,)] 

kqp 

x [f(t k ) - f(t k + eV)} 6{eV + t k - t p _ q + uj q ) (2) 

where e is the electron charge and V is the bias [9]. The differential conductance of 
the junction with antiparallel magnetization directions is just proportional (for zero 
temperature) to the density of the non-quasiparticle states, 

dl/dV oc N NQP (eV) 
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Note that the value N^Qp(eV) vanishes for \eV\ < hu>o, where fiujQ is an anisotropy 
gap in the magnon spectrum [9], which is small, but could be changed by suitable 
substitution [16]. It is worthwhile to note also that the nonquasiparticle states play an 
important role (together with the classical effect of the smearing of the gap by spin- 
disorder scattering) in the depolarization of the states near the Fermi energy at finite 
temperatures; again, estimations of the tunneling magnetoresistance based on a simple 
Stoner-like picture of the energy spectrum appear to be completely unrealistic [17]. 

III. NON-QUASIPARTICLE STATES AND ELECTRONIC STRUCTURE 

To perform the calculations we used the newly developed LDA+DMFT scheme 
[18]. This method is based on the so-called Exact Muffin-tin Orbitals (EMTO) scheme 
[19,20] within a screened KKR [21] approach, frozen core together with the local spin 
density approximation (LDA). The correlation effects are treated in the framework of 
dynamical mean field theory (DMFT) [22], with a spin-polarized T-matrix Fluctuation 
Exchange (SPTF) type of DMFT solver [23]. The SPTF approximation is a multi- 
band spin-polarized generalization of the fluctuation exchange approximation (FLEX) 
[12,13], but with a different treatment of particle-hole (PH) and particle-particle (PP) 
channels. The particle-particle (PP) channel is described by a T-matrix approach [24] 
giving a renormalization of the effective interaction. This effective interaction is used 
explicitly in the particle-hole channel. Justifications, further developments and details 
of this scheme can be found in Ref. [23]. The spin-polaron process is described by the 
fluctuation potential matrix W aa (icu) with a — ±, defined in a similar way as in the 
spin-polarized FLEX approximation [13]: 



The essential feature here is that the potential (3) is a complex energy dependent 
matrix in spin pace with off-diagonal elements: 



where U m represents the bare vertex matrix corresponding to the transverse magnetic 
channel, x <7 '^ <T if UJ ) is an effective transverse susceptibility matrix and Xo _ ° (^) is the 
bare transverse susceptibility [13]. %lo are fermionic Matsubara frequencies and (m) cor- 
responds to the magnetic interaction channel [12,13]. The local Green functions as well 
as the electronic self-energies are spin diagonal for collinear magnetic configurations. 
In this approximation the electronic self-energy is calculated in terms of the effective 
interactions in various channels. The particle-particle contribution to the self-energy 
was combined with the Hartree-Fock and the second order contributions [13]. To en- 
sure a physical transparent description in the current implementation the combined 
particle-particle selfenergy is presented by its a Hartree YP- H ^ and Fock £( TF ) types of 
contributions: 




(3) 



W a - a (iu) = U m ( X a '~ (T (iuj) - Xo~ a M)U" 



rn 



(4) 



(5) 
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where the particle-hole contribution Y,( pK) reads: 

45 = E (6) 

34<t' 

Since some part of the correlation effects are included already in the local spin- 
density approximation (LSDA) "double counted" terms should be taken into account. 
To this aim, we start with the LSDA electronic structure and replace Y, a (E) by T, a (E) — 
S CT (0) in all equations of the LDA+DMFT method [25], the energy E being relative to 
the Fermi energy and Ep = 0. It means that we only add dynamical correlation effects 
to the LSDA method. 

In our calculations we considered the standard representation of the zinc-blende 
structure with an fee unit cell containing four atoms: Cr(0,0,0), As(l/4, 1/4, 1/4) and 
two vacant sites E(l/2, 1/2, 1/2) and El (3/4, 3/4, 3/4). We used the experimental lat- 
tice constant of GaAs(5.65A) and InAs(6.06A) for all the calculations, and the equilib- 
rium value for the bulk CrAs (a eq = 5.8A), obtained by FLAPW [4]. To calculate the 
charge density we integrate along a contour on the complex energy plane which extends 
from the bottom of the band up to the Fermi level [20] , using 30 energy points. For the 
Brillouin zone integration we sum up a k-space grid of 89 points in the irreducible part 
of the Brillouin zone. A cutoff of l ma x = 8 for the multi-pole expansion of the charge 
density and a cutoff of l max = 3 for the wave functions was used. The Perdew-Wang 
[26] parameterization of the LDA to the exchange correlation potential was used. 

IV. RESULTS AND DISCUSSIONS 

The LDA+DMFT calculations were carried out for three lattice constants: the 
GaAs, InAs and the "equilibrium" one. The corresponding LDA computational results 
agree with previous ones [1,4,27]. In order to evaluate the average Coulomb interac- 
tion on the Cr atoms and the corresponding exchange interactions we start with the 
constrained LDA method [28]. In our case the constrained LDA calculation indicates 
that the average Coulomb interaction between the Cr 3d electrons is about U = 6.5 eV 
with an exchange interaction energy about J = 0.9 eV. It is important to note that the 
values of the average Coulomb interaction parameter slightly decrease going from the 
GaAs (U = 6.59 eV) to InAs (U = 6.25 eV) lattice constants, see table Tab. I. This is 
in agreement with a naive picture of a less effective screening due to increasing of the 
distances between the atoms. 

The typical insulating screening used in the constraint calculation [28] should be 
replaced by a metallic kind of screening. The metallic screening will lead to a smaller 
value of U. Unfortunately, there are no reliable schemes to calculate U in metals, 
therefore we choose some intermediate value of U — 2 eV and J = 0.9 eV. It is 
important to realize that there are no significant changes in the values of the average 
Coulomb interaction for the studied lattice structures, and that the exchange interaction 
is practically constant. Note that physical results are not very sensitive to the value of 
U, as it was demonstrated by us for NiMnSb [14]. 

The LSDA and LSDA+DMFT calculation for the density of states, DOS, is pre- 
sented in Fig. 2. Depending on the lattice constant the Cr and As atoms loose electrons 
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and this charge is gained by the vacant sites. As a result the Fermi level is moving from 
the right edge of the gap as for the case of GaAs substrate, towards the middle of the 
gap for a InAs substrate. The Cr moments are well localized due to a mechanism being 
similar to that of localization of the magnetic moment on the Mn atom in the Heusler 
NiMnSb [29]. Note that the local Cr spin moment is more than 3/xb, the As atom pos- 
sesses a small induced magnetic moment (of order of — 0.3/xb) coupled antiparallel to 
the Cr one. The results are presented in table II. Calculated, DMFT, Cr magnetic mo- 
ments increased in comparison with the LDA results, due to the localization tendency 
of the Cr— 3d states as a consequence of correlation effects [18]. 

According to our calculations, the system remains half-metallic with a pretty large 
band gap (about 1.8 eV) for all the lattice constants compared with the band gap of 
the NiMnSb which is only 0.75 eV [29]. 

In Fig. 2 the non-quasiparticle states are visible for lattice parameter higher than 
the equilibrium one, with a considerable spectral weight in the case of the InAs sub- 
strate. InAs is widely used as substrate material for the growth of different compounds. 
According to our calculations this situation is very favorable for the experimental in- 
vestigation of the formation of non-quasiparticle states as well as the investigation of 
their nature. 

Comparing with the case of the GaAs lattice parameter, for the InAs substrate, the 
most significant change in the electronic structure is suffered by the As-p states. Having 
a larger lattice constant the Cr atom acquires a slightly larger magnetic moment. In 
the LDA calculations, nevertheless, the magnetic moment per unit cell is integer 3/ib- 
Expanding the lattice constant from GaAs to InAs lattice the Cr states become more 
atomic like and therefore the spin magnetic moment increases. This is reflected equally 
in the charge transfer which is smaller for the InAs lattice parameters. A larger Cr 
moment induces a large spin polarization of the As-p states, compensating the smaller 
p — d hybridization, the total moment remaining at its integer value of 3/ib [27]. 

Note that the essential difference of the many-body electronic structure for the 
lattice constants of GaAs and InAs is completely due to the difference in the position 
of the Fermi energy with respect to the minority-spin band gap whereas the self-energy 
characterizing the correlation effects is not changed too much (Fig. 3). The total 
density of states N (E) is rather sensitive to the difference between the band edge E c 
and the Fermi energy Ep. If this difference is very small (i.e., the system is close to 
the electronic topological transition E c — > Ep) one can use a simple expression for the 
singular contribution to the bare density of states, 5N (E) oc y/E — E c (E > E c ). The 
appearance of the complex self-energy E (E) = Ei (E) — iE 2 (E) changes the singular 
contribution as 

5N (E) oc [ ^Zl (E) + £| (E) + Z 1 {E)} 1 ' 2 (7) 

where Z\ (E) = E — E c — Ei (E) (cf. Ref. [30]). Assuming that the self-energy is 
small in comparison with E — E c one can find for the states in the gap SN (E) oc 
E 2 (E) I \/E c — E (E < E c ) . One can see that the shift of the gap edge changes drasti- 
cally the density of states for the same E 2 (E) . 

In the remaining part of the paper we elaborate more on the possibility of the 
practical use of tunable properties of non-quasiparticle states in CrAs grown on different 
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substrates. For most of the applications room-temperature and the stability of the 
ferromagnetic state are the most important prerequisites. The ferromagnetic CrAs 
material might be grown on III-V semiconductors similarly to the zinc-blende CrSb 
[33]. We observed the presence of the NQP states for CrAs lattice parameters larger 
than 5.8A It was found experimentally that at 300K, around this value of the lattice 
parameter a stable solid solution of Gao.65Ino.35As is formed [31]. So, from the practical 
point of view a 65% of gallium in a Ga x Ini_ x As compound would constitute the ideal 
substrate on which the CrAs half-metal would present tunable properties of its NQPs. 
Having the Gao.65Ino.35As substrate CrAs could be a part of the epitaxial III-V structure 
providing an easy way to integrate with the existing semiconductor technology. 

In addition, to determine possible substrates for growth of layered half-metallic ma- 
terials, electronic structure calculations were carried out for lattice parameters in the 
range 5.60A to 6.03A [32]. According to these calculations [32], growth with minimal 
strain might be accomplished in a half-metallic multilayer system grown on InAs sub- 
strate which would be the best choice to evidence the NQP states, since the Fermi 
energy is situated in this case far enough from the bottom of the conduction band. 

The existence of NQP states is a manifestation of the many-body interaction at finite 
temperatures, magnetoelectronic applications being able to measure their presence. 

A very high sensitivity of the minority-electron density of states near the Fermi 
energy to the lattice constant opens a new interesting opportunity to run the current 
in CrAs-based tunnel junction. Suppose we have an antiparallel orientation of the 
magnetizations in the tunnel junction (such as shown in Fig. lc), then the I — V 
characteristic is determined by the density of the non-quasiparticle states. Thus if we 
will influence the lattice constant (e.g., using a piezoelectric material) one can modify 
the differential conductivity. This makes CrAs a very promising material with tunable 
characteristics which opens new ways for applications in spintronics. 

This work is part of the research programme of the Stichting voor Fundamenteel 
Onderzoek der Materie (FOM), which is financially supported by the Nederlandse Or- 
ganisatie voor Wetenschappelijk Onderzoek (NWO). 
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TABLES 





Cr 
(a.u) 


As 
(a.u) 


E 

(a.u) 


El 

(a.u) 


U 
(eV) 


J 

(eV) 


T>GaAs 
n MT 

n MT 

nlnAs 

K MT 


2.565 
2.633 
2.751 


2.688 
2.759 
2.883 


2.565 
2.633 
2.751 


2.688 
2.759 
2.883 


6.59 
6.46 
6.25 


0.93 
0.93 
0.93 



TABLE I. The muffin-tin radii of the Cr, As and the two types of empty spheres placed 
as vacancies of the ZB structure. The constreined LDA values for the Coulomb and exchange 
interactions are presented. 





Cr 


As 


E 


El 


Total 


T 


U 


J 






(Mb) 


(Mb) 


(Mb) 


(Mb) 


(K) 


(eV) 


(eV) 


..GaAs 
^LDA 


3.191 


-0.270 


-0.009 


0.089 


3.00 








..GaAs 
r'DMFT 


3.224 


-0.267 


-0.023 


0.067 


3.00 


200 


2 


0.9 


eq. 
^LDA 


3.284 


-0.341 


-0.018 


0.076 


3.00 








eq. 
^DMFT 


3.290 


-0.327 


-0.024 


0.068 


3.00 


200 


2 


0.9 


. .InAs 

Vlda 


3.376 


-0.416 


-0.025 


0.066 


3.00 








..InAs 
r'DMFT 


3.430 


-0.433 


-0.033 


0.043 


3.00 


200 


2 


0.9 



TABLE II. Summarizing table containing the results of our calculation. CrAs magnetic 
moments corresponding to the GaAs, InAs and the equlibrium lattice constant a eq . For the 
later one the value a eq = 5.8A was used, Ref. [4]. Parameters of the DMFT calculations are 
presented in the last three columns of the table. 
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FIGURES 
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FIG. 1. The tunneling transport between strongly correlated ferromagnets. The density of 
states in the lower Hubbard band (a) is provided by standard current states for majority-spin 
electrons (above) and by non-quasiparticle states for minority-spin electrons (below), the latter 
contribution being non-zero only above the Fermi energy (occupied states are shadowed). 
However, the tunneling is possible both for parallel (b) and antiparallel (c) magnetization 
directions. 
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FIG. 2. Cr density of states calculated in LSDA (dashed line) and LSDA+DMFT (solid 
line) methods corresponding to a temperature of T = 200-fT, average Coulomb interaction 
parameter U = 2 eV and exchange J = 0.9 eV. The non-quasiparticle states are clearly 
visible for lattice parameters larger than a eq = 5.8A, in the unoccupied part for minority spin 
channel just above the Fermi level, around 0.5eV. 
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FIG. 3. Energy dependences of imaginary parts of the electron self-energy Ti2(E), for 
lattice constants of GaAs (a), equilibrium one (b), and of InAs (c): e g down solid line, tig 
down decorated solid line, e g up dashed line, tig up decorated dashed line. 
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